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A new relativistic formulation of MOND is advanced, involving two metrics as independent degrees 
of freedom: the MOND metric to which alone matter couples, and an auxiliary metric Qfj,^. 
The main idea hinges on the fact that we can form tensors from the difference of the Levi-Civita 
connections of the two metrics, C^^ = FJ^^ — F^^, and these act like gravitational accelerations. 
In the context of MOND we can form dimensionless 'acceleration' scalars, and functions thereof 
(containing only first derivatives) from contractions of Og^C^^. I look at a class of bimetric MOND 
theories governed by the action / = -{IGtvG)-^ J[Pg^''^ R + ag^''^ R-2{ggY/^ f{t{)alM{T / al)]d^x + 

Imig^Lu, Tpi) + iMigtiv, Xi)y with T a scalar quadratic in the C^^, k — (g/g)^''^, Im the matter action, 
and allowing for the existence of twin matter that couples to g^^ alone. Thus, gravity is modified not 
by modifying the 'elasticity' of the space-time in which matter lives, but by the interaction between 
that space-time and the auxiliary one. In particular, I concentrate on the interesting and simple 
. choice T oc g'^" (C^s^C^^ — C^^Cx-,)- This theory introduces only one new constant, Oo; it tends 

fSJ ' simply to general relativity (GR) in the limit oq 0, and to a phenomenologically valid MOND 

theory in the nonrelativistic limit. The theory naturally gives MOND and "dark energy" effects from 
^ , the same term in the action, both controlled by the MOND constant Oq. As regards gravitational 

lensing by nonrelativistic systems-a holy grail for relativistic MOND theories-the theory predicts 
that the same potential that controls massive-particle motion also dictates lensing in the same way 
as in GR: Lensing and massive-particle probing of galactic fields will require the same "halo" of dark 
matter to explain the departure of the present theory from GR. This last result can be modified 
with other choices of T, but lensing is still enhanced and MOND-like, with an effective logarithmic 
O , potential. 
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I. INTRODUCTION 

> 

, From the inception of MOND [1] it has been elcar that the paradigm needs buttressing by a relativistic formulation. 
0^ Indeed, efforts to construct such a formulation started shortly thereafter, with the tensor-scalar version sketched in 
[2]. This was the first in a chain of theories of increasing force, culminating in the advent of the tensor- vector-scalar 
, theory (TeVeS) of Bekenstein [3]. Some landmarks along this track are described in [3-8]; see, in particular, the 
■ reviews in [5, 8]. All these theories involve as independent degrees of freedom an Einstein metric, whose free action is 
the standard Einstein-Hilbert action, with additional scalar and/or vector degrees of freedom, with their own actions. 
These scalar/vector degrees of freedom are used to dress up the Einstein metric into the 'physical' metric to which 
matter couples. TeVeS has a version of the nonrelativistic (NR) theory proposed in [2] as a NR limit. 

Another line of relativistic theories that aim to reproduce MOND phenomenology has been propounded in [9, 10], 
based on the omnipresence of a gravitationally polarizable medium proposed in [11]. 

Here I propound a new class of relativistic formulations for the MOND paradigm in the form of bimetric MOND 
d (BIMOND) theories. These came to light as follows: I have recently described [12] a new class of nonrelativistic, 
bi-potential MOND theories, a subclass of which is governed by a Lagrangian density of the form 

^ = + ~ "o-^KV.^ - V^)Va2]} + p(iv2 - 0), (1) 



leading to the field equations 

V- [Ai*(lV0*l/ao)V0*] =47rGp, = /? - 



" + A^/^ 2^ 



= V • [(1 - a-^M')W(l)*] = AnGr^p + ■ {M'W(f>*), (2) 

with ip* — (j) — (j). I also described in detail the requirements from a, /3, and A4{z) that lead to the required MOND 
and Newtonian limits of these theories. In particular, I discussed at length the interesting case (3 + a = {P — 1 then 
normalizes G to be the Newton constant), which leads to the field equations 



A<j)* = 4nGp, = V • [(1 + A4')V0*] = 47rGp + V • {M'^(t>*), 



(3) 
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with Ai' a function of (V0*/ao)^, such that A4'{z) for z oo ensures the Newtonian Umit, and Ai'{z) « z~^^^ in 
the MOND regime 2; <C 1. This is a particularly tractable MOND theory, as it requires solving only linear differential 
equations, with the inevitable MOND nonlinearity entering only algebraically. In all the NR theories above, matter 
couples only to one of the potentials: the MOND potential cj), while (j> is an auxiliary potential, and in the special case 
of Eq.(3) their difference (p* is exactly the Newtonian potential of the problem. 

These NR MOND theories have inspired the construction of closely analogous relativistic MOND theories with 
two metrics as independent, gravitational degrees of freedom, which I begin to investigate here.^ This new class of 
BIMOND theories involve only Uo as a new constant. They tend to general relativity (GR) in the limit a,, 0, which 
is a desirable trait. And, they tend to a MOND theory compatible with MOND phenomenology in tlioir NR limit. 

These theories, like all other relativistic versions of MOND proposed to date, must, I believe, be only approximate, 
effective theories to be derived from some more fundamental picture that underlies them. This is pointed to by the 
appearance of an a priori unspecified function in all these theories. 

The use of two (or more) metrics to describe gravity has a long history. For example, Rosen [13] considered bimetric 
theories, where the auxiliary metric is forced to be flat. More recently, it was found [14] that ghosts appear in a 
large class of bimetric theories (apparently not including the present BIMOND). More matter-of-principle questions 
regarding bimetric gravities are discussed in [15-18], but these authors confined themselves to metric couplings that 
involve only the metrics, not their derivatives, as in the case of BIMOND. 

In section II, I present the formalism underlying the BIMOND theories; in section III, I consider the NR limit of 
these theories, showing how they lead to NR MOND theories; section IV demonstrates how the theories go to GR in 
the limit Oq 0; section V discusses lensing; section VI discusses cosmology briefly, and section VII is a discussion. 



II. FORMALISM 



The NR theories mentioned above involve two potentials, the MOND potential (p felt by matter, and an auxiliary 
one (j>. They point to relativistic BIMOND theories involving the MOND metric g'^^, to which matter couples, and 
which in the NR limit reduces to 0, and an auxiliary metric 

Working with two metrics enables us to form nontrivial tensors and scalars from the difference in their Levi-Givita 
connections F^^ and F^^, 

involving only flrst derivatives of the metrics, which is not possible with a single metric. This is particularly pertinent 
in the context of MOND, since connections act like gravitational accelerations. So, without introducing new constants 
in the relativistic formulation we can write Lagrangian functions of dimensionless scalars constructed from a~^C^^ 
that enable us to interpolate between the GR limit, Gq ^ 0, and the MOND limit, ^ 00. 

The tensor C^^ is related to covariant derivatives of one metric with the connection of the other (more generally, 
they relate covariant derivatives of tensors with respect to the two connections): 

9^u:x = Qc.C^x + 9c^,C^^ %.;x = -9c.C;^ - Qac^C-;^, (5) 



Clp = ^g^'igcp-.p + 5/3p:a - g0.f1:,) = —g^'igcp;? + g^p-.c. - gcp j, (6) 

where the covariant derivative (; ) is taken with the connection F^^ and (:) with F^^. We can form various scalars out 
of C^^ and the metrics. One scalar that will be of particular use to us is based on the tensor 

with the same index combination that appears in the expression for the Ricci tensor 

R^, = T" -T" +r^,r^ -F^rt, (8) 



In these theories the two metrics are independent degrees of freedom. Theories like Brans-Dickc, TcVcS, etc., are also sometimes 
described as being bimetric because they involve two metrics, but those two metrics are a priori related conformally or disformally via 
other degrees of freedom such as scalars or vectors. 
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and in R^^ constructed similarly from g^^. One finds 

-R^^ - R^„ = C^^.^ - C^^.^ - T^^. (9) 

Thus, using well known manipulations, the scalar T = g^^T^^ connects the two Ricci scalars R = g^"R^^ and the 
mixed R^ = g'"'R^^ by 

R-R^ = -r + g-'/'ig'/'g'-c;:,),. - (10) 
Similarly, interchanging the roles of g^^ and g,^^, 

R-R^ = -f- g-^i\g^i^rci,). + r^i\9'i^rci:),,, (n) 

where T = ^''"T^^, R^ = ^""iJ^^, R is the Ricci scalar of 5^^, and g and g are minus the determinants of g^,^ and g^^ 
respectively. 

We can construct gravitational Lagrangian densities using the scalars i?, R, R^, and scalars constructed by 
contracting powers of C^,^ with the two metrics and their inverses (there arc also gj g. u) = g''"g,,„, etc. that can be 
used). If we only contract with and <?'"', a quadratic scalar is a linear combination (possibly with coefficients 
depending on scalars such as 5/^ or w) of the following scalars 

ff-C;,C\, C-C„ g^.C^C-^, g^-'C^a, g^.g^^g-"'C;Cl^, (12) 

where C'' = g^^C^^, = C"^. The choice of scalars to be used may be forced on us by various theoretical and 
phenomenological desiderata (see below). The main point is that C^^ have only first derivatives of the metrics, that 
they reduce to derivatives of the potential difference in the Newtonian limit (in the sense to be discussed below) , and 
that we can form dimensionless quantities from them with the MOND acceleration (or a MOND length i = c^/tto). 
As regards the four curvature scalars, is will be advantageous to include them in the action only linearly, and eschew 
terms such as in the voguish f{R) theories. Such nonlinear terms render the theory a higher derivative one, which I 
would like to avoid. Another reason to avoid such terms in the MOND context is that they do not naturally lead 
in their NR limit to a single constant Oo controlling the dynamics.^ Neither obstacle appears if we allow functions of 
scalars made of C°. . These contain only first derivatives of the metrics, and give NR limits in which only a„ appears 

(see below). We see from Eqs.(lO) that g^^^R and g^^'^R^ differ by T plus a total derivative so it is enough to include 
one of these in the action, as we anyhow permit functions of T. The same is true of the pair g^/^R^ and g^^^R. 
Because the number of possible combinations it too large to explore here, I limit myself to the subclass of actions of 
the form 

I=-Y^ J W'^R + ag^'^R-2{ggYl^f{K)r^M{e^T^rHd''x + I^{g,^^^^^ (13) 

where I = c^/ao, k = {g/gY^^-, /(I) = 1, and T^™^ are scalars formed by contracting a product of (even) m C^^, which 
can be used in principle. In what follows, I shall confine myself to quadratic scalars.* I have included two matter 
actions: The first, involves the matter dcgrccis of freedom with which we interact directly, designated symbolically 
as ipi. It contains only the MOND metric to which matter is coupled in the standard way. The other, /„, involves 
other matter degrees of freedom, Xi-> and only p^^, to account for the possibility that controls a matter world of 
its own. There are no direct (electromagnetic, etc.) interactions between the if) matter and the twin x matter.^ 



^ For the same reason I avoid scalars that arc higher order in the curvature tensors, such as the different possible contractions of R^^i, 
with itself or with ij^ii,. These are even less appealing as explained in [19]. 

^ For example, to account for dimensions correctly, a function of R has to be introduced as f{£^R), with i some length scale. The NR 
limit of R includes c~'^A(f> as the dominant term in and second order ones such as c~'*(V(^)^, and c~^(pA((>. Thus, in the argument 

of / the second term will give (V0/ao)'^, with ao = (? jt-, which fits well into the MOND frame. But the first, dominant, term would 
involve a time scale ijc, not an acceleration. When R appears linearly, the first term becomes immaterial in the action, as a complete 
derivative, and we are left with terms that are welcome in MOND [the (/>A(/> term is also (Vi/))^ up to a derivative]. 

■* The MOND constant ao is normalized so that the mass-asymptotic- velocity relation is MGao = V*. It defines the scale length £ that is 
used in the coefficiont and the argument of M . Any dimensionless fa<;tors can be absorbed in the definition of M so that its coefficient 
is = c^^Oq and its argument is as prescribed here. 

^ To obviate possible confusion, note that the twin matter is not to play the role of the putative dark matter in galactic systems; this is 
still fully replaced by MOND effects; see below. 
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I make two requirements of the action: a. Require that it gives a NR MOND theory in its NR hmit. This means 
the following: given a non relativistic system of slow masses one can express the metrics solution of the relativistic 
theory in terms of potentials so that the equations of motion for slow particles in the resulting (multi) potential theory 
are those required by NR MOND, with the appropriate MOND and Newtonian limits; this is a phenomenological 
requirement (by itself it does not dictate the effects on massless particles-e.g., gravitational lensing-even in NR 
systems), b. Require that the action gives GR in the limit Co — ?• 0. This is not a phenomenological necessary (for 
example, TeVeS docs not satisfy it), but I feel that it is highly desirable for various reasons. This automatically 
causes the theory to agree with the stringent constraints from the solar system and binary pulsars-which are known 
to agree with GR-because the accelerations in these systems arc many orders of magnitude larger than a„. I also 
require this limit lest we have to introduce additional constant (s) to the theory, which has to give GR in some limit 
of its parameters. 

When the two metrics arc conformally related, which might be the case in certain circumstances, g^,, = e^^'^^g^^, we 
have C;, = 2t}_^, = -g"-"!},,, T^, = (l/2)(c/^,5°'^i?,„t?,^ - ^,^^,^), T = (3/2)5^''t?,^i?,,. If we a priori constrain our 
metrics to be conformally related (i.e. vary the action only over such pairs) we get the Brans-Dicke theory with the 
choice Ai{z) (X z (and appropriate choice of the constants and /(k), and possibly using T in the argument of A4). 
With a more general form of Ai{z), we then get the relativistic MOND theory sketched in [2]. 

Without the interaction M term, the theory separates into two disjoint copies of GR. It is important to note that 
as a combined structure, the theory then enjoys a larger symmetry involving separate coordinate transformations in 
the two separate actions. This double symmetry has to be brought to bear when solving the field equations of the 
theory, which now satisfy two sets of Bianchi identities. So, eight gauge conditions can, and have to, be employed. 
It is the interaction that breaks this larger symmetry, as, generically, it is only invariant to application of the same 
coordinate transformation to the two metrics. However, under certain circumstances the interaction is symmetric 
under a more extended set of coordinate transformation, and we must be careful then to employ the larger gauge 
freedom. The above mentioned complete decoupling is an example that, as we shall see in section IV, applies in the 
formal limit Oq — )• of the theory (leading, as we want, to GR). It may also happen, in principle, that the interaction 
term vanishes only in some limited regions of space-time; for example, if the extreme GR limit applies in some regions. 
In this case we must allow for gauge freedom involving coordinate transformations that coincide only outside these 
regions, but not inside them. We shall see another example in section III, where the NR limit of the theory has such 
a partial double gauge freedom. 

A. Concrete simple example 

I shall hereafter concentrate on a simple special case of the class. Some generalizations will be mentioned briefly 
below, in this section, and in section VI. 

In the flrst place, I take to be a function of only one scalar, quadratic in the C^^. In particular, I flnd the scalar 
T defined above a natural choice for this argument, as it has the same structure as the first-derivative part of the 
Ricci curvature scalar (not itself a scalar) 

r<^> ^ g-'{Tl,rt^ - Tl^Tl). (14) 

It is well known that one can replace R in the Einstcin-Hilbcrt action by F'^' and still get GR. Here we can also do 
this, replacing also R by the corresponding F*^', and making M a function of T, which is constructed in the same way 
from C^^ . We shall also see that with this choice of scalar argument the NR limit of the theory is especially simple. 

As a further simplification I take a + /3 = 0. This will yield a particularly interesting and simple subclass of theories, 
which turn out to have the theory (3) as their NR limit for slowly moving masses in a double Minkowski background. 
I then take /3 = 1 for G to be Newton's constant. 

Work in units in which c = 1, and use Oq = to highlight the connection with MOND. Also, anticipating the 
expression for NR limit of T, I take the argument of A4 to be — T/2ag. The relativistic action I then consider is 

/ = -Y^ / b'/'^-5'/'^-2(55)'/V('^)«oA^(-T/2a^)]rf^ar + /M(ff,.,V'0-^M(5..,xO- (15) 

[Using Eq.(ll) we can replace the first two terms by {g^^'^g''" — g^^^g'^")Rf^^ + g^^^T.] I take a mixed volume element 
for the interaction term, with / normalized such that /(I) = 1. Note the change of sign in the definition of the twin 
matter action to match the negative sign for the Hilbert-Einstein action of 
Varying over g'"' and over g"" we get, respectively 



(16) 
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G^^+S^^ = -SttG%^, (17) 
where and are the Einstein tensors of the two metrics, 

= Ru,^ — -^Rg^^, G^^ = Ri^^ — -Rg^^, (18) 

and 77j„ arc the matter energy-momentum tensors (EMT); e.g., SI,^j = —(1/2) J g^^'^T,,„Sg^'' , and S^^ are the 
functional derivatives (one with an opposite sign) of the interaction term with respect to the two metrics: 

-2{ggy/^f{K)alM{-r/2a'',)d^x = J {g^/Hg-^S,. - g^lHg"" S,.)d^x. (19) 

For the present choice of the scalar argument of , we have 

S,^ = K_M'T^^ + {n.M'S^),, - A^g,^, (20) 

5,. = {k+M'S'J.., - Kg,., (21) 

= CI - SIC., + \{C' - C^)g,., (22) 

= llC^,^ + g'''C;,J.,^ - g^^C;XJ-^>^ - <llC.> + \g:J'-C^ - \c^g,., (23) 

A. = -^[«/(«)]'a?M, k^ = ~\n-^f{K)\alM. (24) 



Here, 



c^ = 5"^a, k± = k±7(k), Q: = g"'g.., g*=g.,g'"g,., (25) 



and {fjb...v) = {ii...u + z/.../x}/2 signifies symmetrization over the two indices, 
The tensor 

may be viewed as the EMT of the phantom dark matter (DM) ; whereas the A„ term may roughly be viewed as "dark 
energy". Note that the last term in <S^_^, which contributes \[n_M.'{C^ — C^)].^g^,. may also contribute to the dark 
energy due to its form. Define in analogy with T).^ 

The Einstein tensors satisfy the usual Bianchi identities G" = G" = 0,^ derivable from the invariance of the 
Einstein-Hilbert actions to coordinate transformations. In addition, we have here, for the general action (13), a set 
of four identities following from the fact that the mixed term is a scalar; these read 

S;,. - i^-^Sl^ = 0. (28) 

Given that the matter EMTs are divergence free (for matter degrees of freedom satisfying their own equations of 
motion): 7/^_^ = 7^"^^ = 0, the above identities imply four diff^erential identities satisfied by our 20 field equations. If 
we write these equations as = and Q^,. = 0, respectively, then the four relations 

Ql:. - = (29) 



^ For each tensor indices are raise with the corresponding metric; so, e.g., = g^^Gfta, Gj^ = g'"'Gfja 
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hold identically, and, as usual, deprive us of four equations to account for the fact that the solution can be determined 
only up to a coordinate transformation. This seems to leave us with a tractable Cauchy problem, although this require 
more careful checking/ 

Of course, for solutions of the field equations we do have separately 

s;.,. = s;..^ = 0, (30) 

which can be used as useful constraints of the solutions (only one set is independent). 
Note the useful identities 

a = ^r'Sc^:. = -^5"%.. = 2/^,./« C' = -K-\K^g''-).., = -C' - g'r.^, (31) 

<5:.;A = R.. - \Rr.g,. - (i?.. - \r9,.) - (T,. - \Tg,^), (C' - C'),, = R-R^ + T. (32) 

Identities (32) follow from Eqs.(9)(10). Similar manipulations are possible for <S^^, and the field equation (17). 
Contracting Eq.(16) with ^i^" gives 

R - [K.M'iC' - C'')U - K_M'T + 4A„ = SttGT. (33) 

Contracting Eq.(17) with g'"' gives 

R - [K+M'{^ut'C, -C'- g"'q;C-JU + 4A„ = SirGT, (34) 
where ui = g'"'g^,.- We can thus replace Eq.(16) by 

R,. + K_M'{T,^ - ixg,.) - [«_M'(<5(^^a) - C^U + A„g,. = -87rG(r,. - ^Tg,.), (35) 

and similarly for Eq.(17). We can also use identiticis (9-11) to write these equations in different forms. Equations 
(33) (34) can be used to write possibly useful integral (virial) relations by integrating them over space-time, each with 
its own volume element. 

It was deduced in [14] that under certain assumptions about the theory, bimetric theories generically posses ghosts. 
One of their assumptions was that to lowest order in departure from double Minkowski the theory is a sum of 

Pauli-Fierz actions for the different metrics, which are quadratic in the metric departures. This, however, leads to a 
linear theory in this limit, which is at odds with MOND: MOND phenomenology dictates that at (/^^ = g^^ = rj^^ 
any BIMOND theory (or any relativistic MOND theory for that matter) is not even analytic in the squares of the 
departures g^„ — rj^^, g^^ — r]^„ (where the argument of Ai' in the above version of the theory vanishes, and Ai' 
diverges). It thus remains to be seen if obstacles similar to these are at all relevant to BIMOND, and if they are to 
what extent they are deleterious. 

For conformally related metrics = e'"'"^^^„, we have <S^^ = e^^^^S^^. 



B. Generalizations 



Some generalizations of the above simple theory include the following. 

1. Instead of using T as the argument of A^, we can use other scalars, or several scalar variables. A quadratic 
scalar variable can be written, most generally, as 

s = ofrc"C (36) 



As a result of identities (28) and the Bianchi identities, the four expressions G° + — k^^(G'° + S'") contain only up to first time 
derivatives of the metric, and cannot be used to propagate the problem in time. Instead, the initial conditions have to satisfy the four 
equations Q° — k^"^Q" = 0, and the remaining sixteen field equations, with the aid of four gauge conditions, propagate us in time, and 
insure that these four are always satisfied. 
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where Q'^Y" it buih from g,^^,, their inverses, (5°, and scalars such as k and Co. In this case the A terms take a 
more general form, and so do terms that are second order in the C^^. The only terms in S^^ and S^^ that survive in 
the NR limit, which we treat below, are those involving S^^ and S^^. For these we now have for example 

^::. = ^U^:9^,.-U;'g,M'^, Ur = -2QirC;^, (37) 

which I shall need in what follows. 

For example, taking as the argument of M, —C^^x^t-y/'^'^o instead of —T/2a1, would leave us with only the first 
term in expression (22) for S'^^ , and with the first three terms in expression (23) for S^^. 

2. One can consider more general a, /3 values. 

3. We can increase the symmetry with respect to the two metrics by taking interaction terms of the form M{TT), 
M{T)Mif), etc.. 

4. One can make A4 a function of scalars such as k and w. 
Additional generalizations will be mentioned in section VI. 

III. NONRELATIVISTIC LIMIT 

Consider now the NR limit of the theory derived from the action (15). This limit applies to systems where all 
quantities with the dimensions of velocities, such as v, etc., are much smaller than the speed of light. In the 
context of GR this limit is attained by formally taking c — )• oo everywhere in the relativistic theory. In the context of 
MOND one has to be more specific, since system attributes with the dimensions of acceleration, such as w^/i?, V(/), 
etc., cannot be assumed very small in the limiting process, even though they have velocities in the numerator. We 
want to consider systems, such as galaxies, in which these are finite compared with the MOND acceleration, which 
is also a relevant parameter. The NR limit in MOND is thus formally attained by taking everywhere c — )• oo, but at 
the same time £ ^ oo, so that Uq = (? ji remains finite. 

Take a system of quasistatic (nonrelativistically moving) masses, so that to a satisfactory approximation we can. 
as usual, neglect all components of the matter EMT except %o = P- I also neglect here the possible efi'ects of the 
presence of twin matter.^ First, I consider the system in a double Minkowski background. This is aesthetically the 
most appealing option, which I shall assume. It relics on the possibility that on cosmological scales the two metrics are, 
somehow, maintained the same from some symmetry. There are indeed versions of BIMOND [made more symmetric 
in the two metrics than our simple action (15) is] that have cosmological solutions with = 51^^, either at all times, 
or as vacuum solutions, which might be appropriate for today (see section VI). In this case we have C^^ = for the 
cosmological background, and finite C^^ values occur only due to local inhomogeneities. We can then take locally, on 
scales much smaller then cosmological ones, a double Minkowski background. Departures from this assumption will 
be discussed below. 

Write, then, the metrics as slightly perturbed from Minkowski. Because the source system is time-reversal symmetric 
in the approximation we treat it (neglecting motions in the source), we are looking for a solution for which the mixed 
space-time elements of the two metrics vanish.^ We can then write most generally 

9f,. = rit.<^-'2(l>5^<^ + h^,.-, = r?^. - 2(^J^, + A^,, (38) 

where /lo^ = /i^o = ^^o^ = /i^o = 0- We denote the differences 

5*. = -20*'5m- + (39) 

with (j)* = (j) — (j), /i*^ = hi^^ — hu,^. We wish to solve the field equations to first order in the potentials </>, 0, /li^j/iy 
(Roman letters are used for space indices). 

Note that there is a subtlety here (as in all metric MOND theories) due to the fact that the NR MOND potential 
for an isolated mass diverges logarithmically at infinity; so, strictly speaking we cannot formulate a first-order theory 



This is justified if this matter is nonexistent, or of it is smoothly distributed so its local contribution is negligible, or if there does not 
happen to exist a twin body in the near vicinity of the i/) body under study. 

We do not have to assume this a priori; if we do not, the equations themselves will tell us that there is a choice of gauge in which the 
solution satisfies this ansatz; see the end of this subsection. The ansatz simplifies the presentation, and is justified a posteriori by our 
showing below that such a solution exists. 
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for such an isolated mass assuming (f> ^ I at all radii. However, we are, in any event, dealing with an effective theory 
to be understood in the context of the universe at large, and in this context there are no isolated masses, our approach 
is meant to work only well within the distance from the central mass to the next comparable mass, where we can 
assume the first-order theory to be a good approximation. 

To the required order the only nonvanishing components of C^^ are^'^ 

Qo = ^oi = ^io = ~29oo,i = > 

CU = li9*„.+9k, - 9kJ = Ik.. + K.., - h%,) + - <t^*J,, - <t^*J,,. (40) 

These reflect the same relations between the separate connections with their respective potentials. 
The only nonvanishing components of the Ricci tensors (shown here for are 

i?oo = -A<^, i?,, = ijJ,, - A4>5,„ (41) 

with 

H,, = Ah,, + h_,j - 2/i,(,,,),, , (42) 
where h is the trace of hi, . The nonvanising components of the Einstein tensor are 

Goo = -2A</, + ^H, G„ = ^{H,, - ^HS,,) (43) 

{H is the trace of Hij). The same expressions exist for the hatted and for the starred quantities. We are now 
ready to use these expressions in the field equations (16) (17). We neglect the small cosmological-constant terms (in 
line with our assuming background Minkowski metrics), and note that terms such as T^^ are of second order in 
the potentials, so they can be neglected. Also, S^^ and S^^ are linear in components of the tensor C^^, which are 
first order in the potentials; so everywhere else in these expressions we can take the metrics as Minkowski, so that 
/(k) « 1, Kj. « 1, ^ S^, etc.. Also, for the same reason, the covariant derivatives can be replaced by normal 

derivatives. All in all we get that the two terms involving S^^ and S^^ are equal. Thus, taking the difference of the 
two field equations we get 

G;, = -SwGToo, G* = 0. (44) 

Substituting from Eq.(43) we get 

A<1>* - ^H* = 4nGp, H*. - ^H*5,, = 0. (45) 

Taking the trace of the second part we get H* = 0, and substituting in the first we get 

A(p* = A-KGp. (46) 

We impose for (f)* the boundary condition at infinity (f)* — 0, which establishes it as the Newtonian potential of the 
problem. But the second equation (45) does not, in itself, determine /i*., because G*., like G,,, satisfy three Bianchi 
identities G* ^ = 0, which are the reductions of identities (29) to our case (the fourth identity is automatically satisfied 
for our choice of vanishing mixed elements of the metrics) . 

We have only used one of the field equations (or rather their difference). Now consider the first field equation alone 
in the form (35). Again, neglect the second order T terms, etc. to get 

i?,. + [M'{Sl^ - \s^7i,.)h = -^TrGpS,^, (47) 



1° Because the metric derivatives, connections, and curvature components are already first order, all the metrics that are used to contract 
them can be taken as rj^^. 
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where S'^^ is the NR hmit of S^^, and S' its (four) tracx^ The (Oz) components of the equations hold identically (to 
first order) since iS'^. = 0.^^ The (00) and {ij) components give, respectively, 

-A<p + [M'{S^, + ^5-=)],, = -4nGp. (48) 



-i/,, - A(A<5,, + [M'{S^^ - -5*<5.,)].. = -47rGp6^,. (49) 

Multiply Eq.(48) by 6^^ and subtract from Eq.(49) to get 

Ih,, + [M'{S^,-S':MU = 0, (50) 

which I use instead of Eq.(49). Equation (50) does not satisfy any more identities and thus gives six independent 
equations, which together with the above four unused independent equations, and the remaining freedom to choose 
three gauge conditions, should determine the remaining 13 potentials (f), h^j, h*.. 

It is beneficial to employ three of these six equations encapsuled in Eq.(50) by taking its divergence, taking the 
Bianchi identities for into account, to get 

{M'S:,),,.,=Q. (51) 

These three equations, together with the three independent equations in the second of (45) now involve only the six 
h*. (so we managed to decouple these from cj) and hi,; (p*, which also appears in these equations is already known), 
and can be solved for these. Once this is done (imposing boundary conditions at infinity) cf) is determined from 
Eq.(48) by solving a Poisson equation, and hij are likewise determined from Eq.(50) with the aid of gauge conditions. 
The remaining gauge freedom is associated with coordinate transformations that preserve our assumed form of the 
metrics; i.e., near-Minkowski and stationary (time independent and lacking mixed elements). These are of the general 
form 

t = t', a;' = a;'" + e'(x'), (52) 
with e'(x') first order in the potentials. They do not affect g*^ and so leave h*. and 0* intact, changing only by 

Everything so far is valid for an arbitrary choice of quadratic scalar argument. I now specialize to my preferred 
choice of scalar argument — T/2. In this case we have 

Si, = + i(/i* - r), 5-, = i(/i* + hl_^ - h*,J + i[2(/i* - K^JS,, - h*J^, - h*J4. (53) 

What is special about this case is that the space components <Sy depend only on /i*., not on (j)* . This greatly simplifies 
the solution of Eqs.(45)(51), which has to be h*. = (with boundary conditions /i*. at infinity). 

The fact that h*. = causes the C^^, as given in Eq.(40), to be linear combinations of derivatives of <p*, and the 

argument of M.' becomes a function of (V^*)^. Now Eq.(48) reads 

A(f) = A-kGp + V • {M'[{W(I)* /a^f]W(l>*). (54) 

Equation (50) becomes 

i?.. = 0, (55) 



To sec this note that from Eq.(37) we have §1-^^ oc (2(5J^J'' — Qf]'"')C0^- Now, the NR hmit of the Q tensor is eonstructed only from 
Tji^i, and r)^". This means that its only nonvanishiiig components must have three pairs of equal indices. This means, in turn, that the 
only contributions to Sgj come from C^^ with one or three time indices; but these all vanish. 
These are coupled nonlinear equations, since h*^ appear also in the argument os A4' . 

Since g*^ is already first order, in our approximation the transformation affects it only to zeroth order; i.e., not at all. 

All the above holds if we use both T and T* as variables, because they degenerate into one in the NR limit. Also, note that C° = 0, 

and C' = h*j ,j — (l/2)/i* ; so, adding the scalar g^^C'^C" should not change this conclusion. 
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and can be used with three gauge conditions, to determine hi^-as would be done in GR, where this equation is always 
satisfied. This implies, with the appropriate boundary conditions, that there is a gauge in which h^j = 0, and we work 
in this gauge. 

The matter action for a system of slowly moving masses is, to our present approximation, 

lM-\j p{^'^-(l))d^xdt; (56) 

So the motion of such particles is governed by the potential 0, which is thus identified as the MOND potential. It 
is determined from the field equation Eq.(54) (with 0* being the Newtonian potential of the system), which is the 
QUMOND formulation described by Eq.(3), and discussed at length in [12]. Note that it requires solving only linear 
differential equations. To have the required Newtonian limit we have to have M'{z) ^ for ^; oo (i.e., ^ 0). 
In the MOND regime z ^ 1 we have to have M'{z) « z~'^/'^ to get space-time scale invariance, which is the defining 
tenet of the NR MOND limit [20] (the normalization is absorbed in the definition of ao). 
To recapitulate, we end up with the simple result in the chosen gauge: 

9^. = - 2(1)5^,, %^ = T]^, - 24>5^,, (57) 

with <f>* = (f) — (f) being the Newtonian potential, and the MOND potential (p determined from the QUMOND Eq.(3). 
The relation between the first-order MOND metric and the MOND potential is thus exactly the same as that between 
the first-order GR metric and the Newtonian potential. 

Suppose we have not assumed a priori that the mixed elements of the first-order metrics vanish. This does not 
change expressions (40) for C^^, but we now have additional nonvanishing elements 

C'o, = (/^o^,, - KJf^, CI = -{K,, + KJ/2. (58) 

The (Oi) component of the difference Ricci tensor is R*^ = {Ah*^ — /i*^. . ^.)/2. So Eq.(45) is now complemented by 

{AK, - K^,J = 0. (59) 

With our boundary conditions h*. — > at infinity the solution of this equation is h*. ~ v*. for some w*(x) (the left 
hand side is identically divergence free, which leaves us with only two independent equations). Note now that the first- 
order limit of the theory, discussed here, enjoys a less obvious symmetry beside the general invariance to simultaneous 
coordinate transformations: it is invariant under a 'small' transformation of the form t = t' — u(x') applied separately 
to the and the sectors (with u first order in the potentials). In other words, there is symmetry to transforming 
got — 9oi + u^i, Qoi — >■ goi + u,i (hence g*^ — )■ gl^ -I- m* ; u* = u — u), with u and u free for us to choose (the EMTs are 
unchanged to lowest order). We thus have the freedom to choose a gauge for g^^ alone, in which v* = 0, and hence 
/i*. = 0. This means that C^. = C°. = 0, and so the (Oi) components of Eq.(47) read 

(A/io.-/io,„.,) = 0, (60) 

whose solution is hoi = v^i for some w(x). We still have the gauge freedom to chose v{x) = 0, and so we do. We thus 
end up with a gauge in which the field equations themselves dictate hot = hoi = 0, as we assumed a priori. 

The double gauge symmetry we use can be seen to apply directly to the first-order equations. It can be traced 
back to the fact that the NR limit of T is invariant to it: If we do not assume a priori that the mixed elements /i*. 
vanish, then the addition to the lowest order expression for T is C^^Cg^, because C^. is antisymmetric in i,j, while 
C°. is symmetric. But, under the double gauge transformation /i*. changes by w* , so C*^ is invariant, and so is T. 

It remains to be checked if this symmetry is a remnant of some symmetry enjoyed by the relativistic theory itself. 

Anticipating the discussion of the next subsection, note that not all scalar arguments are invariant to this double 
gauge in their NR limit. For example, the change induced in the scalar g^^C'^C'" is — /i*. .Am*. For the more general 
case, Eq.(59) is still valid, and again gives h*^ = v*, while instead of Eq.(60) we have more generally 

(A/io, - ho„„.)/2 + [M'S^,],, = 0. (61) 

The first term is identically divergence free so we can write one of these three equations as 

[M'SI],,,.=0. (62) 

Now, Sg. is linear in v* (which may also appear in the argument of M'), so this equation generically dictates v* = 0. 
For the scalar T this does not work because at this stage we already have S^- = 0, but in return we have the double 
gauge freedom to help us remove v* . For the scalar g^^C'^C" we do not have the double gauge freedom but we can 
write oc 6^Av*, so Eq.(62) gives A{M.'Av*) = 0, which implies v* = 0. In any event we can always have /i*. = 0, 
and continue from there as before to show that there is always a gauge where the mixed elements of the metrics 
vanish. 
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A. Other choices of the scalar argument of M. 

Here I consider the NR limit of theories with other choices of the quadratic scalar argument of M.. The main 
purpose is to sec whether these give theories that are different in their NR limits, and so can be distinguished using 
observations of NR systems such as rotation curves and lensing in galaxies. 

Take then the general quadratic argument as given by Eq.(36). All of our procedures in the present section, up to 
Eq.(52) remain valid. Up to that point we had not made use of the particular expressions for <S^^ and <S^^, only of 
the fact that they become equal to first order in the potentials, and this is still the case.^^ 

Departure from the above occurs, however, for more general scalars in the employment of Eq.(51). Now, the space 
components 5'^. do, in general, depend on the gradient of and it is easy to see that h*. = is no longer a solution, 
in general: Substituting h*- = in Eq.(51) would result in three constraints on the Newtonian potential (p* , which 
it does not satisfy {(j)* is really arbitrary if we allow arbitrary density distributions, including negative ones). So, in 
general, after (p* is calculated as the Newtonian potential of the system, we have to solve the nine coupled equations 
(45) (second part) and (51) (only six of which are independent) for the h*-. After this is done we determine (j) from 
Eq.(48), and hij from Eq.(50) with the aid of gauge conditions. Note that in the Newtonian limit, — )• 0, ^ becomes 
the Newtonian potential, and h^^ — > as fast as M' does. 

We can derive some scaling properties of the h*-, even in the general case: It is easy to see that the second Eq.(45) 
and Eq.(51) are invariant under rrii Xrui, r A^/^r, where nii stand for the masses in the system. This is because 
the only quantities appearing in these equations, beside the variables /i* , are (j)*. and ao both of the same dimensions 
as rriiG/r^. This means that the h*. ,^ are invariant under this scaling. 

Consider now in more detail a spherically symmetric problem. We are looking for solutions in which the various 
potential tensors /i^ , hij , h*. are of the form exemplified by 

h*,=K{r)5,,+K{r)n,n,, (63) 

where n = r/r. It can be shown that h*. satisfying the second Eq.(45), namely, annuling iJ*, is tantamount to 
h* = rhl' , which means, in turn, that h*. = q^ij for some q{r). It thus remains to determine q from Eq.(51). Because 
of the spherical symmetry, the different i components of Eq.(51) give equivalent equations;^^ so we are left with only 
one equation from which to determine q{r). 

Once q{r) is known, we use Eq.(50) to solve for hij. Write in the form (63). We can use the remaining gauge 
freedom to eliminate one of the two functions. In the spherically symmetric case the remaining freedom is to transform 
r = r'[l + e(r')] for some e(r'), treated to first order. This transformation takes h^j — > + 2eSij + 2re'ninj; so, we 
can use such a transformation to eliminate either of the functions in the expression for hij. For example, let us choose 
the gauge in which 

K,=ip{r)6,,. (64) 
The general NR MOND metric is thus diagonal for this choice of gauge with 

ffoo = -l-2</., =<5,,[l-2(</. + ^)]. (65) 

For the form (64) of hij we have 

H,, - iff^y = -(^" + r- V')^.. + (</'" - r- V')"i". = aS,, + bn,n,. (66) 
Note that a' = —r~'^{r'^b)', from the fact that the expression is divergence free. We now use 

F., - \h5^, = -2iM'S:^),,, (67) 

obtained from Eq.(50), to solve for ip. Since the right hand side of this equation is already known to be divergence 
free, from Eq.(51), we get only one independent equation of the form 

r{r^^ip'y = p{r), (68) 



This follows from the asymmetry of to interchange of the matrices, and from the fact that to first order we can put everywhere else 

This equation has to read in the spherical case P[q(r)]r = where P is a differential operator acting on q{r), and we get one equation 

P[q(r)] = 0. 
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where p{r) is a known function. This we finaUy solve for ip, which wc permit to behave asymptotically as ln{r). 

Here we note already an interesting difference from the theories that have T as scalar argument, which have Eqs.(2) 
and (3) as their NR limit (even with general a, /3-see the next subsection). In such theories, in the spherical case the 
MOND acceleration is an algebraic function of the Newtonian one, with the relation being unique for the theory. In 
the general case this is not so: To get the MOND acceleration in the spherical case we apply the Gauss theorem to 
Eq.(48). The expression we then get is some functional of q{r) that cannot be written as a function of the Newtonian 
acceleration —dcj)* /dr. This can lead to different predictions even for massive-particle motions. 

The spherical problem can be solved analytically for the case where (jf = ^r^-for example when we are outside 
the mass where 9 = —1, and we arc in a region where Af'(z) cx for example, in the dcep-MOND regime where 
we will have to have a = 1/4. Then the solution can be shown to be of the form q = Xr^(j)*, with A determined from 
Eq.(51) depending on and a. With this ansatz M' oc (|A|/ao)"^''(a + bX + cA^)-'^r-^'^(^-^), and Eq.(51) then gives 
{A/ao)-^''{a + bX + c\'^y^A{a + bX)r^n = 0, with C = (1 - 2cr)(6l - 1) - 2 (= -3 for the above examples), and a, b 
depending on 6, a and the choice of scalar argument of M for the specific theory (a comes from the terms in 5^ linear 
in the gradient of <p*, and b from those linear in the gradient of /i*.). So A = —a/6 gives us the solution (for the choice 
of T as argument d = 0). Equation (68) then gives tp = ^{A/ao)~^'^ Ar'-'^^ for ( ^ —3, and tp = ^{A/ao)~'^'^ Aln{r) for 
C = —3, with the dimensionless ^ determined. 

Take, for instance the interesting case where we are asymptotically outside matter and in the deep-MOND regime, 
where A = —MG and C = —3. We then get ip = ^{MGaoy/^ln{r), asymptotically. In this case we also have 
(j) = {MGaa)^/'^ln{r), by definition; so, ip = ^0. For example, for the choice of argument —g''''C^^C^^/2ag, I find, 
following the above procedure for the dcep-MOND (ct = 1/4), asymptotic {6 = —1) case: ^ = 4/3. 

In summary, the asymptotic form of the MOND metric is diagonal, with 

5oo = -1 - 20, g,, = S,,[l - 2<P{1 + 0]. (69) 

Remember that while (p — £^(j> in the MOND regime, in the Newtonian regime becomes the Newtonian potential 
while (fi vanishes as fast as dictated by the vanishing of A^' at high values of its argument. 

B. General a j3 values 

The NR limit of the field equation in a theory governed by the action (13) for general a and (3 values is 

PiR.. - lv,.R) + (M'5;J.. = -87rGpS^,S^o, (70) 



a(4. - Iv^.Ii) - {M'S;J,^ = 0, (71) 

with all quantities taken to first order in the potentials (p, cj), h^j, h^j. S^^^^ is linear in the first-order expression for 
C^^, as given in Eq.(40), and M' is a function of a quadratic scalar built from these. Multiply the second equation 
by — /3/a and add to the first to get 

PiK - h"-^*^ + ^i^'K^).' = -87rGpS,oS.o. (72) 

The (Oj) components of this equation hold identically, as before. Equations (41)-(43) are then used to write the (00) 
and (ij) components of this equation as 

- iff*) - ^{M'S'J, = 4nGp, H* - \h*5,, + ^i^t^(A^'5^.).. = 0- (73) 

Taking the trace of the second and substituting for H* in the first we get 

/3A</)* - ^^[M'{Si, + SIJ], = A-KGp. (74) 

Equation (73) comprises seven independent equations that can be solved for (j)* and the /i*.. Once these are known 
we can get cf) and /ly from Eq.(70), or equivalently from 

+ [M'{Sl^ - \s^rj,^)h = -'inGpS,^. (75) 
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Its (00) component gives 

A</, = A-KGp-^p + p-^[\M'{Si, + SU)], (76) 

from which (p can be gotten (since the right hand side is now known). The {ij) component gives [after making use of 
Eq.(76)] 

H,, - \h5^, = ~2r\M'S':X.- (77) 

Prom these hij can be gotten. Note that the divergence of Eq.(77) is identically the same as that of the second 
Eq.(73); so we only have here three independent equations for the six hij. However, again we have the gauge freedom 
to eliminate this indeterminacy. 

Specializing to our preferred choice of scalar S = — T/2 the NR limit of the theory again greatly simplifies since 
from Eq.(53) 5^ does not contain the derivatives of (p* . This implies that the second term in the second Eq.(73) is 
linear in the derivatives of h*., and so the solution of this equation is easy to get: h*. = 0, (again, with the asymptotic 

boundary conditions h*. — >■ 0). This means that the argument oi M' is now a function of 4>* / ao)"^ , and that, in fact, 
»Sj^ =0. As a result the first of Eq.(73) becomes identical with the first of Eq.(2), while Eq.(76) becomes identical 
with the second of Eq.(2). In addition, from Eq.(77) we get /i,, = as before. 

We thus end up with a NR limit in which = r/^^ — 20(5^„ and = ry^„ — 2(/)5^„ as in Eq.(57), with the MOND 
potential ^ determined now from the NR MOND theory described by Eq.(2). This theory has been discussed at 
length in [12]. The relation between the first-order MOND metric and the MOND potential is thus, again, the same 
as that between the first-order GR metric and the Newtonian potential. 

Note in this context as well that using a scalar argument that is a combination of T and T = gi^^C^C", leads to 
the same first-order metric and NR limit. 



C. Other backgrounds 

So far I assumed that the two metrics have the same cosmological background and so, for systems small on the 
cosmological scale both can be taken as nearly Minkowski. 

Here I consider some possible departures from this assumption. One possibility, for example, is that for today's 

cosmology the two metrics are conformally related = Ag^„ with constant A; this still gives C^^ = for the 
cosmological background. In this case we can take locally the background metrics to be g^^ = r/^„, g^^ = Ar?^„ and 
expand around these: 

g^. = v^.-'^(l>^^. + h^,,, g^,. = X{v^. -24>S^^ + h^^), (78) 

instead of Eq.(38). The expressions of F"^, C^^, and i?^,^ in terms of the potentials remain the same as before. Also, 
'small' coordinate transformations of the type shown in Eq.(52) still do not affect the potential differences only the 
hij. Work with the general theory for arbitrary a, /3. The NR limit of the field equations is now 

P{R,. - Iv^.R) + Xf{X-'){M'S;J, = -8nGZ., (79) 

a(4. - lv,.R) - f{X-'){M'Sl),, = 0, (80) 

with the A powers coming from factors such as k±, etc.. Defining a = Xa, and Ad' = Xf{X~^)M' , we get back 
the A = 1 case, but with a replaced by a, and Ai' by A4'. As before, with our favorite choice of scalar argument 
S = -T/2a we have hij = in the appropriate gauge; so we get the relation between the MOND metric and the 
MOND potential as before. 

For a more general background we can write the background metrics,locally for a small system, 

5/.- = 9^- = Xin^u - uS^oS.o)- (81) 

This leads to more complex NR limits, which I do not discuss here. 

We see then that the NR limit of the theory as applied to system that are small on cosmic scale depends on the 
background metrics. If the relation between the background metrics vary with cosmological time (I assume that it 
does not-see section VI) the application of BIMOND to local inhomogeneities also varies with cosmic time. But I will 
not discuss this possibility further here. 
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IV. GENERAL RELATIVITY LIMIT 

I showed in [12] that the NR theories with (3 = 1, but arbitrary a, have a Newtonian Hmit if M'{z) ^ for 2 — >■ c». 
This carries over to the rclativistic theories. The same property of Ai causes the rclativistic theory to go to GR in 
the same hmit — )■ 0, because M' — )■ implies T^^,T^^ — >• 0, and since in this case M.{z)/z must also vanish in the 
limit we also have A„, A„ — )• 0. This gives GR, with g^^^ satisfying the Einstein equation with the standard matter 
EMT as source. In this limit g^^ satisfies its own Einstein equation. It decouples from anyway, so it affects matter 
neither directly nor indirectly. 

We do not have tight phenomenological constraints on how fast MOND approaches Newtonian dynamics for high 
accelerations. There are indications from solar system constraints [1, 21, 22] that M.'[z) decreases at least as fast 
as z^^, but it might turn out to do so much more precipitously. If so, any departure from GR might be practically 
wiped out in a very-high-acceleration system, such as a laboratory on earth, the inner solar system, or a close binary 
pulsar.i'^'^® 

In the limit Oq — >■ 0, the vanishing of T^^ and might occur much faster than that of A„ and A„. For example, 
if A4(oo) is finite, the latter vanish as a^, wile the former might vanish much faster. We can thus, as an intermediate 
approximation, keep the A„ and A^ terms in the theory, and write the limiting field equations (16)(17) as 

R^.- - \R9^.' - Am(oo)ff^. = -SttGT;., (82) 

- \r%^ - A„(oo)5^„ = -SwGZ^. (83) 

For general P values it remains to be checked whether the requirement on M' from the NR limit, deduced in [12] 
suffices to guarantee a GR limit for — )• 0. 

V. GRAVITATIONAL LENSING 

One of the two main phenomenological duties we expect from a rclativistic MOND theory is to predict gravitational 
lensing correctly. In particular, we know that lensing analysis of galactic systems indicate mass discrepancies that are 
not very different to those derived from massive-particle motions (e.g., rotation curves). In other words, lensing in the 
MOND regime is found observationally to be greatly enhanced over the GR prediction without DM. Reproducing this 
fact has been a pressing desideratum in constructing relativistic MOND theories, achieved finally in TeVeS through 
the efforts of Sanders [4] and Bekenstein [3]. 

In the present BIMOND class such enhanced, MOND-like lensing is predicted naturally by all the theories in the 
class. For a choice of the scalar argument that is a combinations of T, T*, and T (and with any a, /3), relation 
(57) between the first-order MOND metric and the MOND potential holds. So, to this order the MOND connection 
r^_^ is expressed in terms of the MOND potential (j) in the same way as the GR connection is expressed in terms of 
the Newtonian potential. This means, in turn, that the MOND potential describes the dynamics of both massive 
and massless particles in the same way as the Newtonian potential does in GR. In other words, such theories predict 
that analyzing lensing and massive-particle dynamics by a NR system assuming GR, should give the same effective 
potential (or the same distribution of "phantom matter"). This is consistent with observations. 

We also saw that there are other choices of the scalar argument of A4 for which the NR MOND metric is characterized 
by additional potentials hi^. However, these vanish quickly for accelerations much above a„, while in the MOND regime 
they are of the same order as the MOND potential (p. We then expect these theories to yield somewhat different lensing 
to that expected with the GR relation between metric and potential, albeit still with the MOND characteristics. 

For example, we saw that far from a central mass M, and in the deep-MOND regime, the form of the MOND 
metric is given by Eq.(69), and this leads to lensing that is multiplied by a factor 1-1-^/2 over that expected from the 
MOND potential with the GR prescription. 

Eventually, by comparing lensing and massive-particle dynamics in the low-acceleration fields of galaxies or other 
galactic systems, we may be able to differentiate observationally between theories with different scalar arguments. 



The presence of the galactic field still induces the departures discussed in [21]. 

A theory like TcVcS also has a surrogate for ao that appears in its NR limit, which is constructed out of the constants characterizing 
the theory. However, TeVeS does not become GR in the limit ao — !• and this leaves possibly detected effects even in an isolated solax 
system, the binary pulsar, etc., even with the very high accelerations characterizing them [3, 8, 23—25] 
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VI. COSMOLOGY 

I cannot at present offer a specific BIMOND cosmology. There are two obstacles to doing this. In the first 
place wc cannot even pinpoint the exact BIMOND theory out of the various versions possible. NR phenomenology 
can assist somewhat in pinpointing the NR limit. However, even for a given NR MOND theory there are different 
relativistic versions having this limit, which differ greatly in the relativistic regime, and specifically in their application 
to cosmology. 

Second, as we well know from a century of experience with GR, having the underlying theory is one thing; pinpointing 
the cosmology is another: In dealing with standard GR cosmology, the cosmological evolution and the present state of 
the universe does not emerge uniquely from first principles. There are major observational constraints, assumptions 
about symmetries, initial conditions, and matter content, that are put in by hand into cosmological theory. We do 
not know, for instance, the initial conditions for our universe from first principles, so an initial singularity (as opposed 
say to a steady state universe with continuous matter creation, or to a static universe as Einstein would have it 
initially) is imposed by hand. Early infiation (the mechanism for which is still moot) is put in by hand to account 
for various observational facts. Cosmic acceleration, whose cause remains unknown, is imposed by hand by invoking 
'dark energy', modified gravity, or other mechanisms. The material content of the universe (e.g. the very existence of 
baryon asymmetry-mechanism still unknown) is an input in cosmology. A priori, we could have had a cosmos with a 
space that is inhomogeneous or anisotropic on cosmological scales; but, the cosmological principle is imposed based 
on what our eyes tell us about our universe. All this is even more acute in light of recent developments in quantum 
cosmology. 

In the case of BIMOND we are on even shakier ground when coming to construct a cosmology. Here we are dealing 
with two space-times, only one of which we can sense directly. We have no direct knowledge of many of the global 
properties of the other space-time. Did it have a Big Bang? Did it undergo inflation? Is it spatially flat? (in Rosen's 
theory the auxiliary metric is constrained to be flat)? Is it spatially homogeneous and isotropic? What is the nature 
of the twin matter? Is it there at all? Is it always homogeneously distributed, or does it clump? Is it characterized 
by the same baryon asymmetry, etc.? 

MOND phenomenology in systems that are small on cosmological scales is particularly simple and clear cut in a 
double Minkowski background. Such a background applies if in cosmology = 5^^. BIMOND cosmology that have 
disparate metrics, so that the cosmological values of C^^ are appreciable, do not seem to make phenomenological 
sense in small systems. 

I thus assume, as an additional cosmological assumption to the many above, that on cosmological scales we have 
= This could emerge, as an external constraint on BIMOND, from the world picture that underlies it. Or, 
more appealingly, it could, at least, correspond to a solution of BIMOND itself in some version. 

I thus consider here briefly only cosmologies with = g,^,^, or with the somewhat relaxed assumption = Ag,,,^, 
with a constant A. In either case we have C^^ = in cosmology; so, finite values of C^^ are produced only due to local 
inhomogeneities. This greatly simplifies the equations of motion, since the only contributions of the interaction term 
that survive are the A„ terms. With this ansatz, the equations of motion for the more general action (13) are then 

aG^, + qalM{0)g^. = -8nGZ^{Xg^^, Xi), (84) 

where 

q = {X/2)[Kf{K)]U.., q = -(A-V2)[/c-V(/^)]U-i, (85) 

and I used the fact that with the above ansatz G^^ = For the two equalities in equation (84) to hold simultane- 
ously, we need to start with a BIMOND theory with some symmetry with respect to the two metrics. As an example 
consider a special, more symmetric, case of the gravitational Lagrangian in Eq.(13) written as 

Pg^/^{R - 2alM) + ag^/^{R - 2alM). (86) 

This choice corresponds to /(k) = {pK + aii~^)/{a + P), and M. = {a + j3)M.. It gives q = P/{a + 13), q = Xa/{a + p). 
So the field equations are now 

+ alM{0)g,^ = -STrGl3-'Z.{g,u,^i), 



G^^ + XalM{0)g^^ = -8TrGa-%^iXg^^,Xi)- 



(87) 
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Clearly, = (A = 1) always corresponds to a vacuum solution of this theory, with both space-times being a de 
Sitter or Anti de Sitter, with a cosmological constant A = —a-^M{0). Furthermore, if we also have from symmetry, 
for two identical configurations in the two sectors, = {a/ 13)%^,^^ then the two equations are the same even 
with matter. The cosmology we then get is, quite interestingly, standard GR cosmology (taking (3 = 1) with A as 
cosmological constant. This would be reassuring since it would automatically insure that we are not bereft of the 
successes of standard cosmology regarding inflation, nucleosynthesis, etc.^° 

This picture would also force us to consider more seriously the nature of the twin matter, and its possible visible 
effects in our space-time. If it is homogeneously distributed it will be difficult to detect any direct effects of it. If 
it clumps it could have various effects; for example, it may produce some effects that are otherwise attributed to 
cosmological dark matter. 

It is also possible that BIMOND can replace cosmological DM by the distribution and fluctuations in T^^, which 
is constructed from the two metrics alone, not directly from matter. This would be similar in vein to what has been 
discussed in connection with such a possible role of auxiliary fields in other theories [8, 18, 23, 26]. 

The above is only one example of a BIMOND theory that has a cosmology with g^^ = This particular version 
should not be assumed for the special case (3 + a = 0, since then Ai drops from the theory in the NR limit. However 
there are other versions of BIMOND that can accommodate our cosmological ansatz for this case. For example, we 
can take as the gravitational Lagrangian density (say with /3 = 1) 

g'/\R - 2pal) - a''\R - 2pal) - 2{ggf'^ f{K)alM = 

= g'/^R - g'/^R - 2a'',{ggnK - pn-^ + f{K)M], (88) 

This theory also has the NR limit we discussed above, and has a cosmological solution with g^,^ = g^^ if a certain 
relation between p, p, and Ai (0) holds. For example, if A^(0) = we have the symmetric theory with p = p, in which 
case —pal is the cosmological constant. More generally, we can take the gravitational Lagrangian density (keeping 
the symmetry) 

g'/'{R - 2Mal) - g'/^R - 2Mal) - 2{ggf'^ f{K)alM = 

= g'/^R - g'/^R - 2al{ggnK - k-')M + finjM] (89) 

(with M. also a function of the scalar argument T/a^). Then, M does not appear in the NR limit on a double 
Minkowski background, which is governed by A4, and a cosmology with our ansatz is a solution, with —a'lM.{0) being 
the cosmological constant [for Ai (0) = 0] . 

We saw that in the GR limit a'^M^ plays the role of a cosmological constant [M^ = M{oo)], while in the present 
context it is a^MiO). More generally, the A„ oc a^M term, which is variable, and possibly other terms in T^^, may 
give rise to "dark energy" effects. Ai{z) need not change much for the full range of z, since for small values it also 
has to go to some constant, as Ai{z) « Mq + (4/3)2:"^/^ for z <C 1. In fact, changes in Ai{z) over the full z range are, 
generically, of order unity, since M.^^ — A4{z) = J~ M.'{z)dz is of order unity if M.' decreases beyond z = 1 faster 
than z~^ (unlike A^, which is known only up to an additive constant, M.' is determined by MOND phenomenology). 
Aloo is a dimensionless constant characterizing the theory. If \Ma^ \ ~ 1, then |A^| is always of order unity. We then 
automatically get the well known, but otherwise mysterious, proximity between a^, as determined from the dynamics 
of small systems, and A, the density of 'dark energy', as deduced from cosmology.^^ 

Note that irrespective of the relevance to MOND, bimetric theories of the type presented here provide a frame for 
discussing 'dark energy' as modified gravity, which could be an alternative to schemes such as f{R) theories. 



This ensures that without the interaction, physics is the same in the two sectors. 
This hue of thinking seems to indicate that the theories with a = /3 = 1 are preferable. 

The possible proximity of A and af^ , thus hinges on the dimensionless M. being of order unity. Because of the way the normalization 
of M is defined, this means that the scale over which M varies as a function of T, and the scale that determines the magnitude of the 
M term in the Lagrangian, which have the same dimensions of length"^, are also of the same magnitude. This need not be the case, 
just as not all mass parameters that appear in the standard model of particle physics are of similar values. So, the apparent proximity 
A ~ Og that we get here is only a plausibility not a corollary. 
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A. Deep-MOND relativistic systems? 



In principle, BIMOND theories enable one to study the structure of deep-MOND, relativistic systems such as 
dccp-MOND black holes. As has been stressed many times in the past, such a deep-MOND system would have to 
have its typical curvature radius much larger than the MOND length f = (? ja^- This length is, however, of the 
order of the Hubble radius today, and certainly in the past. In practice then, the universe seems to be the only such 
low-acceleration (rather, intermediate-acceleration) relativistic system, at present. 



VII. DISCUSSION 



I have described a class of bimetric MOND (BIMOND) theories. Matter lives in the space-time described by one 
of the metrics, which, in turn, couples to another through the interaction M. term. If we heuristically view gravity 
as reflecting an effective 'elasticity' of space-time we can view the double-metric nature of our theory as representing 
two coexisting elastic bodies, each with its own elasticity as encapsulcd in the respective _R, R terms in the action. 
Thus MOND departure from GR is introduced not through a modification of the 'elasticity' properties of space-time, 
but rather through the interaction of the space-time that is the arena for matter with the auxiliary one. The strength 
of the interaction between these two space-time 'membranes' depends on the gradient difference. The response of 
our home space-time to matter is affected by its interaction with the other space-time, which modifies its effective 
elasticity. However, once its shape is determined, this home space-time affects matter in the standard way. With 
our assumption that on cosmological scales = the two membranes are, in a sense, stuck together on these 
scales, and 'separate' only locally due to inhomogeneities. Such heuristics may help pinpoint the fundamental concept 
underpinning the MOND paradigm. For example, it may give meaning to the length I = (? ja^ that appears in the 
NR limit as a^. 

The BIMOND theories have the (yet unproven) potential to account for all the components of the dark sector 
(galactic DM, cosmological DM, and dark energy) from one term in the action, all controlled by a^. 

My main objective has been to point out that there exists such a class of relativistic theories that have MOND-like 
theories as their NR limit, and, which produce enhanced, MOND-like gravitational lensing. We are, however, still far 
from pinpointing the exact version of the theory that is the most suitable. This is particularly true in the context of 
cosmology, which depends crucially on the choice of version. Hopefully, theoretical and phenomenological constraints 
will be brought to bear on this by future studies. It remains to be seen whether a version of BIMOND can be found 
that pass muster given all such requirements. Recent discussions of matter-of-principle questions, such as the causal 
structure of bimetric theories of a different type (where the interaction term is a hmction of the metrics themselves, 
not their derivatives) can be found, e.g., in [16, 17]. 

Finally, it has to be realized that however useful such theories may turn out to be, they must be only effective, 
approximate theories, as evinced by the appearance of the a priori unspecified function M. and the length I (or the 
MOND acceleration Oq). These will, hopefully, be calculated from a theory at a deeper stratum. 
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